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Abstract 

We study neutralino sectors in extensions of the MSSM that dynamically generate the u-term. 
The extra neutralino states are superpartners of the Higgs singlets and/or additional gauge bosons. 
The extended models may have distinct lightest neutralino properties which can have important 
influences on their phenomenology. We consider constraints on the lightest neutralino from LEP, 
Tevatron, and {g — 2)« measurements and the relic density of the dark matter. The lightest 
neutralino can be extremely light and/or dominated by its singlino component which does not 
couple directly to SM particles except Higgs doublets. 
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I. INTRODUCTION 



One of the most important challenges to the Standard Model (SM) of particle physics 
is the cosmological observation that the dominant component of matter in the universe is 
unexplained. The preferred interpretation of the observational data is that this matter is 
massive (cold), electrically neutral (dark), and also stable or long-lived. The Supersymmetric 
extension of the SM provides a well-motivated candidate for cold dark matter. Supersym- 
metry (SUSY) predicts the superpartners of the neutral gauge bosons and the neutral Higgs 
bosons with masses of electroweak (EW) scale and couplings of EW strength. The lightest 
neutralino is the lightest Supersymmetric particle (LSP) over most of the parameter space. 
Assuming i?-parity conservation that prevents rapid proton decay, the lightest neutralino is 
stable. 

The Minimal version of the Supersymmetric SM (MSSM) assumes .R-parity conservation 
and has an extra Higgs doublet and general SUSY breaking soft terms [l[ . Extensive studies 
of the MSSM show that its lightest neutralino has the right ranges of mass and interaction 
strength to be a good cold dark matter (CDM) candidate 0. Although the MSSM may 
be the optimal low energy Supersymmetric model with minimal extension of the fields and 
symmetry, the model may not fully describe the TeV scale physics. It is important to see if 
other versions of the Supersymmetric SM can also give acceptable CDM. 

A good starting point for alternatives to the MSSM is its theoretical weakness. The 
MSSM has a fine-tuning problem, the so called /i-problem The //-term is the only 
dimensionful parameter in the SUSY conserving sector of the MSSM, and it is required to 
have the same scale as the SUSY breaking parameters (EW/TeV-scale soft terms) to give 
an EW scale Higgs vacuum expectation value (VEV) without fine-tuning. The MSSM by 
itself does not explain why the \x parameter should be at the EW scale. Considering that 
the motivation of the supersymmetrization of the SM was to resolve a fine-tuning problem 
(gauge hierarchy problem), the fine-tuning problem of the MSSM is a serious issue. 

We consider various extended MSSM models that resolve this /i-problem and compare 
their lightest neutralino properties to that of the MSSM. In these beyond-MSSM models, 
at least one Higgs singlet is commonly present and generates an effective /i parameter when 
the associated symmetry is broken at the EW or TeV scale. Because of the superpartners 
of the Higgs singlets or the extra gauge boson, the neutralino sector in these models may be 
significantly different from that of the MSSM. We investigate the mass and the coupling of 
the lightest neutralino of each model allowed by the model parameters and the experimental 
data. 

In Section |Hj we describe the models and their neutralino sectors. In Section IIII1 we 
analyze the lightest neutralinos allowed by the model parameters and direct experimental 
constraints. In Section IIV1 we discuss indirect observational constraints on these lightest 
neutralinos. In Section |VJ we discuss our numerical results, and then summarize our results 
in Section IVT1 

II. MODELS 

The extended MSSM models that we consider are the Next-to-Minimal Supersymmet- 
ric SM (NMSSM) Q, the Minimal Non-minimal Supersymmetric SM (MNSSM) a.k.a. the 
nearly Minimal Supersymmetric SM (nMSSM) [5|, the C/(l)'-extended Minimal Supersym- 
metric SM (UMSSM) p, and the {7(l)'-extended Supersymmetric SM with a secluded U(l)'- 
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TABLE I: Higgses and Neutralinos of the MSSM and its extensions 



Model 


Symmetry 


Higgses (CP even 


CP odd, charf 


?ed) 


Neutralinos 


MSSM 












NMSSM 


Z 3 


+ HIAI 






+ s 


nMSSM 




+ HlA» 






+ s 


UMSSM 


u(iy 


+ H$ 






+ s,z' 


S-model 


u(iy 


i fjO uO ttO ttO 


AO AO AO 
A 2 ' ^3 1 ^4 




+ S,Z', Si,S2, s 3 



breaking sector (S-model) |2|. All of these extended models prevent the /i-term (fiHiH 2 ) and 
allow an effective //-term (SHiH 2 ) through a VEV (S) of a Higgs singlet associated with a 
new symmetry. The NMSSM and the nMSSM adopt discrete symmetries while the UMSSM 
and the S-model use an Abelian gauge symmetry spontaneously broken by the Higgs singlet. 

A. Superpotentials 

The superpotential of Higgses (both isospin doublets and singlets) for each model is given 
below 1 . 

Hmssm = ^HiH 2 (1) 
HWssm = h s SH l H 2 + —S H^mssm = h s SHiH 2 + aS (2) 

H^umssm = h s SHiH 2 H^s-modei = h s SH\H 2 + \ S S\S 2 S3 (3) 

The other parts of the superpotentials are the Yukawa terms of the MSSM and possible extra 
terms related to the exotic chiral fields in the U{1)' models needed to cancel anomalies. The 
exotic chiral terms are model-dependent and we do not specify them here, assuming that 
the masses of exotic fields are heavy enough to give insignificant effects to the EW scale 
phenomenology that we are interested in. Specific examples of models with exotic fields can 
be found in Ref. 

The Higgses and neutralinos of each model are listed in Table H] In the £/(l)'-extended 
model, the addition of one Higgs singlet does not give an additional CP odd Higgs since a 
goldstone boson is absorbed to be the longitudinal mode of the massive U{1)' gauge boson, 
Z'. 

We use the common notation of h s for the coefficient of the SH\H 2 term in each model 
for easy comparison. In every model the dynamically generated effective /1 parameter is 
given by 

/"off = K (S) (4) 

and therefore the VEV of the Higgs singlet or the symmetry breaking scale needs to be at 
the EW/TeV scale. 



The term aS in the nMSSM is a loop-generated tadpole term that breaks the discrete symmetry; see Ref. 

a. 
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B. Neutralino mass matrices 



The MSSM has 4 neutralinos (B, W 3 , H%, H%) while the NMSSM has 5 neutralinos 
(MSSM components + S). The neutralino mass matrix of the NMSSM in the {B, W3, Hf, 
S} basis is given by 



n 2i 



( 



\ 



Mi 


-9iVi/2 
g x v 2 /2 






M 2 
g 2 v x /2 
-92V2/2 




-giVi/2 

g 2 vi/2 


— A*eff 
-HeftV 2 / S 



9iV 2 /2 
~92V 2 /2 



-fleSVi/s 



\ 



-HcSV 2 / S 
-HeSVl/s 

2ks/V2 ) 



(5) 



The first 4x4 submatrix corresponds to the MSSM limit, which can be obtained from Eq. 
(JSJ by taking s ^> 0(EW). Depending on the value of the parameter k, the 5th component 
may be very heavy and decoupled from the other EW scale components or almost massless. 
The NMSSM assumes a discrete symmetry Z 3 to avoid the /i-term, but allows S 3 term in the 
superpotential. The VEVs of the Higgses are defined as (Hf) = y= with \Jv\+v 2 = ■ 
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GeV and (S) = -A^. The gauge couplings are g\ = e/ cosOw and g 2 = ej smOw- The NMSSM 
is one of the simplest extensions of the MSSM, but the Z3 symmetry predicts domain walls 
which are not observed 

The nMSSM p| was devised to avoid the domain wall problem while maintaining a 
discrete symmetry. The nMSSM has the same 5 neutralinos and the same mass matrix as 
Eq. (jSJ) except that the k term in the (5, 5) entry vanishes since it is not allowed by the 
discrete symmetry. 

The UMSSM uses an Abelian gauge symmetry instead of a discrete symmetry and thus 
is free from the domain wall problem. With the superpartner of the U(l)' gauge boson (Z'), 
the UMSSM has 6 neutralinos (NMSSM or nMSSM components + Z') and its mass matrix 
in the basis of {B, W 3 , H%, H$, S, Z'} is given by 



/ 



\ 





-g\vi/2 

giv 2 /2 







M 2 
g 2 Vi/2 
-g2V 2 /2 







-givi/2 

g 2 v x /2 


— A'eff 
-He&V 2 /s 



gxv 2 /2 
-5^2/2 

— /^eflf 








-fi eS v 2 /s g Z 'Q' Hl v i 

-fletfVx/s g Z 'Q' H2 v 2 



\ 







9Z'Q' Hl Vl 9Z>Q'h 2 V2 9Z'Q' S S 



gz'Q's 



(6) 



The first 5x5 submatrix corresponds to the nMSSM limit (or NMSSM limit in the special 
case k = 0) that can be obtained from Eq. (jUJ) by taking M\> 3> 0(EW). In this limit 
the mass of the Z'-ino becomes very large and this component decouples from the others. 
Here gz> is the U{1)' gauge coupling constant and Q' is the U{1)' charge. The U{1)' charges 
should satisfy 



Q'h, + Q'm ^ Q'h, + Q'h, + Q's 







(7) 



'H ± ^ ^H 2 T- u ^ ^H 2 

in order to replace the /x-term with the effective /z-term dynamically generated by the Higgs 
singlet S. For the numerical analysis in this paper we use the 77-model charge assignments 
and the Grand Unification Theory (GUT) motivated gauge coupling, g z >: 



Q' Hl 



1 



2V15 



Q'i 



H 2 



2VT5 



9z> 



9\ 



(8) 
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The S-model was introduced to resolve tension between the EW scale /z e g and the heavy 
Z' (up to multi-TeV scale). It is basically the extension of the UMSSM with 3 additional 
Higgs singlets to provide additional contributions to the Z' mass while keeping /j e g = h a (S) 
at the EW scale. The S-model has 9 neutralinos (UMSSM components + Si, S 2 , S3), and 
its mass matrix has 3 more columns/rows added to Eq. (jBJ). (See Ref. 0] for the S-model 
and its full 9x9 neutralino mass matrix.) Its first 6x6 sub matrix corresponds to the 
UMSSM limit, which can be realized by taking $1,2,3 ^ O(EW) with X s comparable to 
gauge couplings. However, the most realistic case is for small X s and large ,§1,2,3 |7|, in which 
case four of the neutralinos, consisting almost entirely of Z', Si, S 2 , and S 3 , essentially 
decouple from the others. Since the full 9x9 matrix has a number of free parameters, we 
consider only this decoupling limit when we discuss the light neutralinos, where there are 5 
neutralinos with masses and compositions the same as the nMSSM 2 . 



C. Interesting limits for the neutralino masses 



In this section, we analyze the neutralino mass matrix more closely, especially for the 
nMSSM and UMSSM limits. The diagonalization of the neutralino mass matrix is accom- 
plished via a unitary matrix N as 

N T M x oN = Diag(M x? , M x o, - • •)■ (9) 

The singlino (S) composition of the lightest neutralino (x?) is |iV 15 | 2 . 

The nMSSM allows the possibility of a very light or massless neutralino which is mainly 
singlino. This is apparent from (0) in the limit of very large s (and k — 0), for which one 
obtains a massless singlino. However, the tendency for a light singlino-dominated neutralino 
persists even for smaller s. In fact, it is possible to have an exactly massless eigenvalue, 
which occurs for 

Det(M x o) = (^p) 2 (M|(Mi cos 2 9 W + M 2 sin 2 9 W ) - ^iM x M 2 sin2/3) = 0, (10) 
which leads to 



Aff « 0.8/iM 2 sin 2/3, (11) 

where we have assumed the gaugino mass unification condition M\ ~ O.5M2. Eq. (jlljl can 
easily be satisfied with EW/TeV scale values of // and M 2 . It is interesting that this can 
occur only for /iM 2 > 0, which is favored by the deviation of the muon anomalous magnetic 
moment data from the SM expectation ^(j 01 , and also by b — > 57 data in most cases 0] • 
In practice, parameters satisfying condition (jllj) are often excluded by the chargino mass 
and Z width constraints discussed in Section UTTl Nevertheless, there are allowed points lying 
nearby in which the combination of moderate s and the smallness of M§ — 0.8/iM2 sin 2/3 



2 The four decoupled neutralinos typically consist of one heavy pair involving the Z' and one linear com- 
bination of Si, S2, S3, as well as two more states associated with the orthogonal combinations of Si, 2, 3 
0,0. The latter can be light, or even be the lightest neutralino in limiting cases. We do not consider 
that possibility here. 
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lead to a very light singlino-dominated neutralino. For large enough s the neutralino can be 
massless. We will further discuss the massless neutralino scenario later. 
Now we consider interesting limits of the UMSSM. 

(i) In the case that only My, the diagonal element of the 6th row/column of Eq. JBJ), is 
very large compared to other mass parameters, the Z' will be very heavy and decoupled, 
leaving a mass matrix similar to that of the nMSSM. 

(ii) When only s is very large compared to other mass parameters (with /i g at the EW 
scale, which requires h s ~ 0), the elements of the 5th row/column in Eq. © are small except 
for the 6th component. Then the effective mass matrix for the 5th and the 6th components 
is 

M ^Lks 92 f 3 )' (i2 » 

which gives almost degenerate physical masses of gz'\Q's\ s - These two heavy states approx- 
imately decouple from the other four neutralinos, which are MSSM-like. 

(iii) When both My and s are large compared to other mass parameters and p e g is at the 
EW scale (i.e., h s is very small), the effective mass matrix for the 5th and 6th components 
is 



which has eigenvalues 



WT±V Mi , and i±vT±V Mi , (M) 

where p = \gz'Q's s /My\. The corresponding neutralino masses in various limits are 

0, \M V \ (p<l) (15) 
p\My\, p\M v \ (p>l) (16) 
0.62|Ml/|, 1.62|Mi/| (p = 1) (17) 

The p < 1 case corresponds to taking the UMSSM to the nMSSM limit (My > C(EW)) 
and then to the MSSM limit (s ^> O(EW)) consecutively. Then a very light neutralino 
state along with a very heavy neutralino state occurs. Besides these lightest (dominated by 
singlino) and heaviest (dominated by Z'-ino) states, the other neutralinos are basically the 
same as the MSSM neutralinos. The p ^> 1 case is similar to the Eq. ()12j) case. For the 
range of p = 10 ~ 0.1, the mass ratio of the lighter state to the heavier one is 1 ~ 100. 

2 

(iv) The case of gaugino mass unification My = Mi = \%M<i ~ 0.5M 2 will be considered 
in the next section. 



III. DIRECT CONSTRAINTS ON THE LIGHTEST NEUTRALINOS 

Here we discuss the allowed mass range of the lightest neutralino after incorporating 
direct constraints from the LEP experiments and the model parameter structure. 
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A. LEP constraints on the light chargino mass and the Z width 



All the models considered have only one charged gaugino (W^) and one charged Higgsino 
(H^) and thus have a common chargino mass matrix which is the same as in the MSSM, 

M *±={ roM M2 a ^ M - Sin/3 V (18) 

X \ v V / 2M W /COS/3 /i e ff J 

The LEP2 data requires the light chargino mass to be M+ > 104 GeV LJ which gives 
lower bounds on M 2 and // e ff for a fixed tan/3. 

For a lightest neutralino with M x o < M z /2, the LEP constraint r^ xp -T| M = (-2.0±2.6) 
MeV |l4| must be satisfied. Since the Z boson does not couple directly to the singlino or Z'- 
ino, the Z-\\~x\ coupling for every model is the same as that of the MSSM. The expression 
for the partial width for Z decay to a pair of the lightest neutralinos is 

r^ M = 9l i 3lm ^ 4,2)2 («! " (2M x ,)f 2 e(M z - 2M,,). (19) 

For a lightest neutralino of very small mass {M x o <c Mz/2) to be allowed 3 , its two dou- 
blet Higgsino compositions should be either almost identical (|A 13 | 2 « |A 14 | 2 ) or negligible 
(|A 13 | 2 , |A 14 | 2 <C 1) by gaugino or singlino dominance. 



B. Lightest neutralino mass range allowed by direct constraints 

We evaluate the bounds on the lightest neutralino mass (M x o) and the singlino compo- 
nent of x? in the MSSM, the NMSSM, the nMSSM, and the UMSSM with gaugino mass 
unification. The nMSSM also represents the NMSSM in the k — > limit, the UMSSM in 
the My ^> Mi limit, and the decoupling limit of the S-model. 

We require that the tree-level masses satisfy the direct LEP limits of M + > 104 GeV, and 

A 1 

T z ^ x o x o < 2.3 MeV (95% C.L.). Bounds from naturalness and perturbativity constraints 

0, Ell [l^ are a l so imposed on the couplings of 0.1 < h s < 0.75 and \/h 2 s + k 2 < 0.75 (for 
the NMSSM) 4 . The LEP2 Higgs mass bound of > 114 GeV does not apply directly to 
the extended MSSM models where the physical Higgses exist as mixtures of doublets and 
singlets fl^ . 

We choose a phase convention in which \i and the VEVs are real and positive, k and 
the gaugino masses can in principle be complex, but we restrict our considerations to real 
values. We mainly consider the favored case of positive gaugino masses, but comment on 
the effects of negative masses. We scan M 2 , \i = 50 ~ 500 GeV with a step size of 1 



3 There are some points involving two or more very light neutralinos, for which the lightest neutralino 
satisfies the Z width constraint but the others do not. This generally occurs for small M2 and fi values 
which are already excluded by the chargino mass constraint. 

4 An exact h s bound (and its source) may be a little different depending on models, but we assume a 
common bound for easy comparison. A lower bound on the NMSSM |k| is fuzzy except that k ^ should 
be satisfied to avoid an unacceptable Peccei-Quinn symmetry. We set |«| > 0.1 as a lower bound; results 
for a smaller \n\ can be described by an interpolation of the nMSSM result (k — > limit). 
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TABLE II: The lower (upper) bounds on the lightest neutralino mass (M^o), the singlino composi- 
tion ( 1 7V"i 5 1 2 ) and the parameter values at the bounds in various Supersymmetric models for several 
fixed tan/3 values. All mass units are in GeV, except for T z _^ o o in MeV. For the UMSSM, the 
77-model charge assignment is assumed. We consider the range 0.1 < h s < 0.75 and only cases 
where the direct experimental constraints M^± > 104 GeV and T z ^ x o x o < 2.3 MeV (95% C.L.) 
are satisfied. The scans are made for M%, fj, = 50 ~ 500 GeV and s = 50 ~ 2000 GeV. For the 
NMSSM, the parameter range \k\ =0.1 ~ 0.75 is scanned, and the k > and k < cases are listed 
separately. 



tan /3 


Model M o AT 15 1 2 
x l 


p M 2 a K 


h 3 M ± Too 
xf z^xlxf 


1 


MSSM 55(242) 
NMSSM 55(242) 0.00(0.00) 
0(242) 0.78(0.00) 
nMSSM 2(83) 0.99(0.71) 
UMSSM 0(242) 0.60(0.00) 


481(500) 121(500) 

481(500) 121(500) 1850(1605) 0.15(0.30) 
242(500) 265(500) 470(1995) -0.10(-0.40) 
142(121) 486(500) 2000(230) 
171(500) 444(500) 365(2000) 


104(420) - ( - ) 
0.37(0.44) 104(420) - ( - ) 
0.73(0.35) 173(420) 0.00( - ) 
0.10(0.74) 124(105) 0.00( - ) 
0.66(0.35) 149(420) 0.55( - ) 


2 


MSSM 55(243) 
NMSSM 48(243) 0.23(0.00) 
17(243) 0.78(0.00) 
nMSSM 1(50) 0.98(0.58) 
UMSSM 26(243) 0.61(0.00) 


486(500) 118(500) 

168(500) 190(500) 320(2000) 0.10(0.65) 
255(500) 500(500) 495(1030) -0.10(~0.20) 
142(151) 248(500) 2000(285) 
228(500) 487(500) 430(1200) 


104(425) - ( - ) 
0.74(0.35) 104(425) - ( - ) 
0.73(0.69) 234(425) 2.29( - ) 
0.10(0.75) 104(136) 0.03( - ) 
0.75(0.59) 208(425) 2.29( - ) 
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MSSM 54(247) 
NMSSM 21(247) 0.43(0.00) 
37(247) 0.71(0.00) 
nMSSM 0(6) 0.98(0.87) 
UMSSM 39(247) 0.55(0.00) 


468(500) 110(500) 

202(500) 133(500) 405(2000) 0.10(0.65) 
216(500) 476(500) 415(980) -0.10(-0.20) 
161(355) 321(500) 2000(670) 
173(500) 499(500) 330(1180) 


104(441) - ( - ) 
0.71(0.35) 104(441) 2.29( - ) 
0.74(0.72) 206(441) 2.30( - ) 
0.11(0.75) 145(335) 0.04(2.26) 
0.74(0.60) 165(441) 2.30( - ) 


50 


MSSM 53(248) 
NMSSM 16(248) 0.48(0.00) 
39(248) 0.66(0.00) 
nMSSM 0(10) 0.98(0.76) 
UMSSM 41(248) 0.50(0.00) 


451(500) 108(500) 

181(500) 133(500) 405(2000) 0.10(0.65) 
191(500) 475(500) 365(980) -0.10(-0.20) 
500(246) 500(118) 2000(590) 
147(500) 500(500) 280(1170) 


104(445) - ( - ) 
0.63(0.35) 104(445) 2.29( - ) 
0.74(0.72) 184(445) 2.30( - ) 
0.35(0.59) 445(104) 0.04(2.29) 
0.74(0.60) 143(445) 2.30( - ) 



GeV 5 , s = 50 ~ 2000 GeV with a step size of 5 GeV, tan/5 = 0.5, 1, 1.5, 2, 10, 50, and 
\k\ = 0.1 ~ 0.75 with a step size of 0.05 (for the NMSSM). The gaugino mass unification 

t- „2 

relation My = M 1 = |^M 2 ~ 0.5M 2 is assumed. The Higgs singlet VEV s cannot be too 

small in the UMSSM because it is responsible for the mass of the Z' , but there are some 
possible ways to get around this, as we will discuss in Section IIVBI 
The bounds that we obtain on the the lightest neutralino mass are 

53 GeV < M x o < 248 GeV [MSSM] (20) 

GeV < M x? < 248 GeV [NMSSM] (21) 

GeV < M x o < 83 GeV [nMSSM, S-model (decoupling limit)] (22) 

GeV < M x o < 248 GeV [UMSSM] (23) 

Table |H] shows the lightest neutralino mass {M x o) and its singlino composition (|V 15 | 2 ) 
as well as other relevant parameter values at the bounds. The results in the table are 
only for positive gaugino masses, but negative gaugino masses do not change these ranges 
significantly 6 . 



5 We exclude very small or ijl to avoid very light non-singlino states. These are also excluded by the 
chargino mass constraints. 

6 With M 2 = -50 500 GeV, the X °i mass ranges are M % o = 39 - 254 GeV (MSSM), M x o =0-254 
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Figured] shows the M x o dependence on tan (3 and | iV ls | 2 (and also |iV"i 6 | 2 for the UMSSM). 
The dashed lines are the MSSM bounds which neglect the rather weak tan/? dependence and 
use the values of Eq. ([20)1 . The mass bound dependence on tan/? is quite sensitive in some 
models. For example, at tan/? « 1, the NMSSM and the UMSSM have massless Xi while 
the nMSSM has an upper M x o bound; the MSSM violates the LEP2 constraint at this 
tan f3. Singlino dominance is typical in the Xi °f the extended MSSM models. Especially, 
when the M x o is much smaller than the MSSM lower limit (M x o ~ 50 GeV), the singlino 
is always the dominant component. We will further discuss the allowed mass ranges along 
with additional indirect constraints in Section |V| 



IV. INDIRECT CONSTRAINTS ON THE LIGHTEST NEUTRALINOS 



A. Cold dark matter relic density and muon anomalous magnetic moment 

The cold dark matter density is tightly constrained by the WMAP (CMB) and SDSS 
(Large Scale Structure) data jl8( to be (with la uncertainty) 

ft CDM /i 2 = 0.12 ±0.01 (WMAP + SDSS) (24) 

where h = 0.72 ± 0.08 is the present day Hubble constant H in units of 100 km s -1 Mpc -1 
[lit ]. We assume that Xi is the sole dark matter and impose the relic density constraint 7 of 
Eq. J2H). 

Another important experimental result that can constrain new physics models is the BNL 



E821 measurement of the muon anomalous magnetic moment = (g — 2)^/2 [111] . The 
deviation from the SM is 



Aa M = cv(exp) - a^(SM) = (23.9 ± 10.0) x 10~ 10 (25) 

when the SM prediction is based on the hadronic contribution from e + e~ data. However, 
the 2.4cx deviation is reduced to 0.9<r if indirect hadronic r decay data are used instead for 
the SM prediction. 

The dominant Supersymmetric contributions to (g— 2) M come from the chargino-sneutrino 
and the neutralino-smuon loops. The Aa^ result practically constrains the sign of M 2 to 
be positive in our sign convention of fi e s > 0. The predicted (g — 2) M value sensitively 
depends on tan j3 and the scalar muon mass. The muon trilinear scalar coupling dependence 
is neglected in our calculation of a M . 

Typically, a large tan/3 value is favored to explain the sizable 2.4<r deviation, but with 
a rather small scalar muon mass, quite small tan/3 values are also acceptable. Figure |2Ja) 
shows how the acceptable parameter ranges from Figure E^a) change due to the Aa^ con- 
straint for two choices of scalar muon mass, ttil = ttie = 100 GeV (dots) and 500 GeV 
(squares). 



GeV (NMSSM), M % o = 0.4 ~ 96 GeV (nMSSM), M x o = 39 ~ 254 GeV (UMSSM). 
7 It was emphasized in |20j that predicted values of f^cDM^ 2 lower than the observed range in (|24|l may be 
allowed if there are nonthermal production mechanisms. 
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FIG. 1: Scatter plots of (a) the Xi mass (M^o) versus tan/3 and (b) the Xi singlino composition 
( | iVi 5 1 2 ) versus M x o for various models. The crosses in (b) represent the Z'-ino composition ( | iV"i6 1 2 ) 
for the UMSSM. The direct constraints of Section ing are imposed. The solid curves represent a fixed 
set of inputs M 2 = 250 GeV, \i = 250 GeV, s = 500 GeV and k = 0.5 (for the NMSSM) without 
the direct constraints. The upper (lower) UMSSM singlino band in (b) corresponds approximately 
to moderate (large) value of s, with the lightest neutralino being MSSM-like for large s. 



B. Z' boson mass 



In U (l)'-extended MSSM models, the mass of the new gauge boson Z' is also an important 
constraint. The Tevatron Run2 dilepton data places a lower bound on the Z' mass of 
500 ~ 800 GeV, with the exact bound depending on the model [2). In the UMSSM, the Z' 
mass is given by 

Ml, = g%, {Q ,2 Hl vl + Q'l % v\ + Q>*s 2 ) . (26) 

Accordingly, the value of s should be at the few TeV level to satisfy the experimental Mz> 
bound; this high value of s, with /i g ~ 0(EW), requires very small h s which implies a 
fine-tuning. 

The bound on s can be significantly reduced if there are additional contributions to Mz>- 
In the S-model, for example, the Z' mass gets additional contributions from 3 more Higgs 
singlet VEVs si,2,3, which practically removes any lower bound on s if the Sj's are at the 
' j_ (^^/"— scQjls 

M% = 9% {QhA + Q'lA + Qs* 2 + EliQsA) (27) 
These multiple singlets help to keep /i e g = h s -^ at the EW scale (resolving the //-problem) 
without fine-tuning even for a multi-TeV scale Z' . It is also possible that Z' could be 
leptophobic, in which case the Tevatron bound on Mz< is greatly reduced. Then s is not 
severely bounded even if there are no additional contributions to Mz>- 
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tan|3 M 2 (GeV) 
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FIG. 2: (a) Similar plot as Figure ^a) that satisfy the 2.4<r allowed {g — 2)« measurement with 
rriL = rriE = 100 GeV (dots) and 500 (GeV) (squares), (b) M x o versus Mi with extended ranges 
of M 2 and fi up to 1000 GeV. The solid curves are for (j, = 250 GeV, s = 500 GeV, tan (3 = 2 and 
K = 0.5. 



V. DISCUSSION OF THE LIGHTEST NEUTRALINO 
A. In the MSSM 

The lightest neutralino mass bound in the MSSM for various conditions has been well 
studied [22]. In our parameter range, the upper bound is 248 GeV, which is associated with 
the assumed upper bound on Mi. The Xi mass bound increases if the Mi bound increases 
(Figure 12(b)). The lower bound on M x o is determined by the light chargino mass bound of 
M x ± > 104 GeV; the bound is M x o > 53 GeV (roughly half of the chargino mass bound), 
for which the T z ^ x o x o constraint is irrelevant. In this M x o range the XiXi annihilation 
cross section is of suitable size for the general MSSM to provide an acceptable relic density 
through multiple annihilation channels j2|. The (g — 2) M deviation can be simultaneously 
explained without difficulty |23|. The LEP2 SM-like Higgs mass bound of > 114 GeV 
excludes tan/5 ~ 1 in the MSSM, though it does not apply to other models with a Higgs 
singlet. 



B. In the NMSSM 

Systematic studies of the NMSSM neutralinos can be found in Ref. |24|. The upper 
bound on the Xi mass is 248 GeV, as in the MSSM. At this point the singlino component 
is absent and k is rather large, 0.65 for k > (—0.20 for k < 0). With a large enough 
us value the singlino component can be very heavy and decoupled, leaving the rest of the 
matrix similar to the MSSM. 
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FIG. 3: (a) Illustration of the massless neutralinos (dark curve) and nearby light neutralinos 
(M o < 1 GeV) (orange region) in the M2-H plane along with constraints from Tz, M ± and h s 
in the nMSSM or the S-model in the decoupling limit. Here tan/3 = 10 and s = 700 GeV are 
assumed. The gap in the Tz exclusion region is due to the emergence of a very light gaugino for 
Mi 2 ~ 0. (b) The relic density for the nMSSM or the S-model with Z-pole annihilation. The solid 
curve is for a fixed set of values of [i = 200 GeV, s = 400 GeV and tan/? = 1.5 with M2 varying 
over ~ 500 GeV. 



With positive k, the lower M x o bound is 16 GeV. With negative k, a massless Xi occurs 
at tan/? = 1 with a singlino composition of 1 7V"i 5 1 2 of 0.78. Since the massless Xi state 
occurs as the result of mass matrix mixing, loop effects such as threshold corrections on the 
gaugino masses do not prevent the appearance of the massless state. Its mixing with the 
right-handed neutrinos may result in interesting neutrino physics (25^ . 

This model is disfavored by the non-observation of cosmological domain walls predicted 
by the discrete symmetries of the model . However, there is an approach to interpret the 
domain wall as the dark energy |26| . The domain wall network has to be strongly frustrated 
(nearly static) to satisfy the CMB isotropy constraint. In that case the equation of state of 
dark energy is expected to be close to w — —2/3, which is disfavored by SNIa data and also 
by joint analysis with other cosmological data |27J but is not excluded. 



C. In the nMSSM or the S-model 

In the nMSSM (also in the S-model in the decoupling limit), a vanishing lower Xi mass 
bound occurs with dominant singlino composition 1 7V"i 5 1 2 ~ 1, where T z ^ x o x o of Eq. (JTHJ) is 
negligible. As discussed after Eq. (|llj). the condition for a massless neutralino can be easily 
satisfied regardless of the tan/3 and s choices. However, the massless condition for small 
tan/3 violates the M x ± > 104 GeV constraint. This is why the lower bound is not exactly 
zero for small tan f3 in Table ITU but can nevertheless be very small for large s. Large s is also 
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s (GeV) s (GeV) 



(a) (b) 

FIG. 4: Scatter plots of (a) M x o and (b) | 1 2 versus s. The solid curves are for M2 = 250 GeV, 
fi = 250 GeV, tan (3 = 2 and k = 0.5. The Mz 1 bound s > 1500 GeV is approximate and can be 
evaded, as discussed in Section HV Bl The crosses in (b) are the Z'-ino composition (|iVi6| 2 ) in the 
UMSSM. 



needed for a light Xi to ensure singlino domination. Figure Efa) illustrates how the massless 
or very light neutralinos of M x o < 1 GeV can appear without violating direct constraints. 

The M x o upper bound is 83 GeV at tan/? = 1, which is considerably lower than those of 
other models 8 ; the bound is due to the upper limit imposed on h s , and the bound significantly 
increases without this constraint. Figure^a) shows that M x o decreases with s, as we expect 
from Eq. (|T0|) . As seen in FigureEfb), the singlino component in Xi increases with increasing 
s and becomes dominant for sufficiently large s. While M x o increases with decreasing s, 
s cannot be too small because of the h s < 0.75 requirement. The singlino composition 
at the maximum M x o is still dominant with | iVi 5 1 2 of 0.7. When there are other sizable 
components in x\i the V z ^ x o x o constraint is not easy to escape unless tan/5 ~ 1 where the 
Z-XrXi coupling (Eq. (fTTIjl ) vanishes. At tan/3 = 1 (or v\ = v 2 ), the contributions of the 
two doublet Higgsinos to Xi are identical (up to sign) as Eq. (jSJ) suggests. 

Because of the small Xi mass in this model, most annihilation channels for the MSSM y? 
relic density calculation become irrelevant, and the Z pole is the dominant channel 
Figure Efb) shows the neutralino relic density in this model through the Z pole annihilation. 
The direct constraints of Section IIHI are applied. To reproduce the acceptable relic density 
with only the Z pole annihilation contribution, the lower M x o bound 9 is M x o ^ 30 GeV, 
while the upper bound remains the same as discussed above. This is the most severe lower 
bound on M x o in the model. However, the Higgs masses are not bounded by the LEP2 data 
because of possible mixing between Higgs doublets and singlets, and very light Higgses may 



8 Including values with < 0, the M % o upper bound can increase, but only up to 100 GeV. 

9 The bound can be lowered a bit for a weaker r 7 __o v o constraint. 

^ ' X l X i 
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provide sufficiently large annihilation so that even lighter neutralinos are allowed by the 
relic density constraint. For an explicit calculation of such a light neutralino relic density 
through a light pseudoscalar state, see Ref. We also refer to Ref. |3(| for interesting 

physics associated with a light pseudoscalar Higgs boson. 

As can be seen in Figure ^a), the lightest neutralinos that are massive enough for the 
annihilation through the Z pole (M x o ^ 30 GeV) are allowed only for small tan (3. This raises 
concern since the Supersymmetric contribution to (g — 2)^ is proportional to tan (3, and a 
large value of tan f3 is favored to explain the considerable 2.4<r deviation of the experimental 
value from the SM prediction. Nonetheless, a common solution was found to exist that can 
explain both the acceptable relic density through the Z pole and the deviation of (g — 2)^ 
in this model, as illustrated in Figure |5Ja) j3l[ . 

D. In the UMSSM 

In the UMSSM, the upper bound on the lightest neutralino mass is typically given by 
the maximal value of Mi (with s ^> Mi), where Xi is essentially the MSSM-like neutralino 
with almost no singlino composition. The Xi mass bound does not increase above 420 GeV 
when the M 12 ranges are extended, unlike the MSSM case (Figure I2^b)). This is because 
of our restriction s < 2000 GeV: for the extended range one is in a regime similar to 
Eq. (|13|). with M x o given by the smaller of the eigenvalues in Eq. (fT4^) . As the solid curves 
of Figure 0] illustrate, both the mass and the singlino composition of the lightest neutralino 
increase with s until s reaches a certain value where that neutralino is no longer the lightest 
and the singlino composition of Xi plunges to zero. 

A vanishing lower bound on M x o occurs at tan (3=1 and s = 365 GeV with a singlino 
composition of | iVis | 2 ~ 0.6 and a Z'-'mo composition of | iVie | 2 ~ 0.2. Though the singlino 
component does not saturate Xi here, the next-to-dominant component, Z'-ino, also does not 
contribute to the Z decay width, and such a light neutralino survives the T z ^ x o x o constraint. 

The Tevatron lower bound on Mz> is too large for the Z' to make a significant contribution 
to (g — 2) /t throu gh a Z'-loop, though the fi-Xr^j contribution is modified by the Z'-mo 
component in x\ |31| . Figure |2fa) shows the parameter space constrained by (g — 2)^ for 
two values of the smuon mass. 

A value of s ~ 0(EW) can be rather problematic, as discussed in Section II V Bl unless 
either there are additional contributions to Mz* (S-model) or the leptonic coupling is small 
(leptophobic Z' model). If we take a TeV scale lower bound on s (for example, in the s ^ 1.5 
TeV range in Figure 0J) to satisfy the Tevatron Mz> bound of 500 ~ 800 GeV, the singlino 
and Z'-mo components are negligible and the lightest neutralino mass becomes similar to 
that of the MSSM. 

VI. SUMMARY AND CONCLUSION 

Although Supersymmetry at the TeV scale is well-motivated, the MSSM is just one of 
the possible realizations. In fact, the theoretical /^-problem suggests that the MSSM is 
incomplete. The solution to the //-problem suggests that an appropriate direction to extend 
the MSSM is to have an extra Higgs singlet whose VEV gives the effective /i-term of the EW 
scale. Extensions of the MSSM have extra neutralinos, and the composition of the lightest 
neutralino involves extra components beyond those of the MSSM. Because of this, both the 
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FIG. 5: (a) An illustration of a common solution of the ^ x oh 2 (dark dots, 3<r of WMAP+SDSS) 
and the {g — 2)^ (cyan region, 95% C.L. of E821) in the M 2 -fi plane along with other constraints 
in the nMSSM or the S-model. The Z-pole annihilation and mi = tue = 100 GeV are used. 
The gap in the Tz exclusion region is where | iV"i3 1 2 ~ 1 7V"i4 1 2 . (b) The allowed mass range of Xi 
after applying direct constraints (model parameters, gaugino mass unification My = M\ ~ O.5M2, 
ilf^i, Ar^) and indirect constraints ((g — 2)^, Q x oh 2 , Mz', domain wall). The M^o bounds in the 
nMSSM are intended to be illustrative and are not necessarily quantitatively precise. 



mass and couplings of the lightest neutralino are modified from the MSSM. The lightest 
neutralino (x?) is interesting both in particle physics (as the LSP) and cosmology (as the 
CDM), and it is therefore important to study and compare properties of the Xi i n extended 
MSSM models. 

We explored Xi properties in various extended MSSM models. We examined constraints 
from the experimental bound on the Z — > X1X1 contribution to Tz, the lower bound on 
M+, and constraints from perturbativity and naturalness of the SHiH 2 and S" 3 coupling 

Al 

strengths. We also considered constraints from the relic density of the CDM (Q x oh 2 ), the 
experimental deviation of (g — 2) of the muon from the SM expectation, and the Tevatron 
lower bound on Mz>- Distinguishing properties of the lightest neutralino arise in extended 
MSSM models, such as distinct x\ mass ranges (Eq. (f2T?|) - (j2HJ)), frequent singlino dominance 
(Figure ^b)), importance of the Z pole annihilation channel in the relic density calculation 
(Figure Efb)) and different tan/3 dependences of the upper and lower mass bounds (Figure 
ma)). Some of the extended MSSM models can be considered as limits of the other models. 
For example, as far as the neutralino sector is concerned, we can consider that MSSM C 
nMSSM C UMSSM C S-model. 

Approximate lightest neutralino mass ranges in the models considered are illustrated in 
Figure Efb). The dashed regions are disfavored by the indirect constraints. After the Mz> 
lower bound is imposed, the UMSSM bound becomes similar to the MSSM. In the case of 
the nMSSM there is a tension between the (g — 2) M constraint which favors small M x o (or 
large tan/3) and the relic density constraint which favors large M x o (or small tan/5). 
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The properties of the CDM particle, even if it is the lightest neutralino, may be quite 
different from the MSSM prediction. For instance, it could be extremely light and/or domi- 
nated by the singlino, which does not directly couple to SM particles except Higgs doublets. 
Similar distinctions of models may occur in the Higgs sector. 

Even if a low energy Supersymmetry is correct, its realization may depend on the model. 
The measurement of the mass of the lightest neutralino and the determination of its couplings 
will be particularly useful in testing the MSSM and its extensions at colliders. 
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